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Titel 

Modellering  bescherming  NBC-kleding 

1  Stromingsprofiel  tussen  kleding  en  huid 

Auteur(s) 

Dr.  ir.  P.  Brasser  en  Ir.  J.  Kaaijk 

Datum 

juni 1999 

Opdrachtnr. 

:  B93KL458 

Rapportnr. 

:  PML 1998-A107 

Er  zijn  diverse  experiraentele  methoden  op  het  TNO  Prins  Maurits  Laboratorium 
(TNO-PML)  beschikbaar  om  de  bescherming  van  NBC-kleding  tegen  de  blootstel- 
ling  aan  cheraische  strijdmiddelen  te  evalueren.  Om  meer  inzicht  te  verkrijgen  in 
de  processen  die  de  beschermende  werking  van  NBC-kleding  bepalen,  wordt,  naast 
de  experimentele  toetsing,  gewerkt  aan  de  modeliering  van  de  bescherming.  Na 
ontwikkeling  van  zo’n  model  kan  dit  worden  toegepast  bij  de  evaluatie  van  nieuwe 
NBC-kleding  en  kan  het  bijdragen  in  de  keuze  van  zo  goed  mogelijk  beschermende 
kleding  met  zo  gering  mogelijke  fysiologische  belasting.  Bovendien  kan  het  wor¬ 
den  gebruikt  bij  de  ontwikkeling  en  verbetering  van  testsystemen  voor  NBC- 
beschermende  kleding. 

Om  tot  een  volledige  modeliering  van  de  bescherming  van  een  compleet  kleding- 
systeem  te  komen,  moeten  de  volgende  processen  beschreven  worden: 

•  het  luchtstromingsprofiel  om,  door  en  onder  de  kleding; 

•  de  penetratie  van  het  chemisch  strijdmiddel  door  het  kledingmateriaal; 

•  het  concentratieprofiel  onder  de  kleding  en  het  damptransport  naar  de  huid. 

In  eerder  onderzoek  is  een  model  ontwikkeld,  dat  de  penetratie  van  toxische  stoffen 
door  NBC-kledingmateriaal  voorspelt.  In  dit  rapport  wordt  een  beschrijving  gege- 
ven  van  een  model  voor  het  luchtstromingsprofiel  om,  door  en  onder  de  kleding. 
Een  lichaamsdeel  wordt  hierbij  gekarakteriseerd  door  een  cilinder,  die  op  een 
bepaalde  afstand  door  kledingmateriaal  is  omspannen.  Voor  het  berekenen  van  de 
luchtsnelheid  door  de  kleding  moet  het  drukprofiel  om  de  kleding  bekend  zijn.  Bij 
het  opstellen  van  dit  model  is  uitgegaan  van  een  drukprofiel  dat  in  de  literatuur 
gegeven  wordt.  De  resultaten  van  het  model  werden  vergeleken  met  die  van  een 
empirische  formule  die  wordt  gebruikt  bij  de  damptesten  van  NBC-beschermende 
kledingmaterialen.  Bij  de  experimentele  condities  waarbij  de  empirische  formule 
werd  bepaald,  bleken  de  resultaten  een  goede  overeenkomst  te  vertonen.  Door 
middel  van  een  parameterstudie  werd  de  invloed  van  diverse  parameters  in  het 
model  onderzocht.  Indien  de  dikte  van  de  luchtlaag  onder  de  kleding  afneemt,  of  de 
diameter  van  de  cilinder  toeneemt,  kan  dit  leiden  tot  een  afname  van  de  luchtsnel¬ 
heid  door  het  kledingmateriaal.  De  grootste  effecten  worden  gevonden  bij  mate- 
rialen  met  een  lage  luchtweerstand. 

Ter  verbetering  van  het  model,  zullen  drukprofielen  om  cilinders  met  kleding 
gemeten  worden  of  numeriek  worden  berekend. 
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Het  model  zal  als  uitgangspunt  dienen  voor  het  modelleren  van  het  concentratie- 
profiel  onder  de  kleding  en  het  damptransport  naar  de  huid.  De  resultaten  hiervan 
zullen  worden  vergeleken  met  de  experimentele  studie  die  in  het  kader  van  het 
project  ‘evaluatie  kledingtesten’  wordt  uitgevoerd  (A97KL423). 


Samenvatting 


Er  is  een  model  ontwikkeld  dat  de  luchtstroming  om,  door  en  onder  NBC- 
beschermende  kleding  berekent.  Een  lichaamsdeel  wordt  hierbij  gekarakteriseerd 
door  een  cilinder,  die  op  een  bepaalde  afstand  door  het  kledingmateriaal  is  om- 
spannen.  Om  het  stromingsprofiel  te  berekenen  moet  het  drukprofiel  om  de  kleding 
bekend  zijn.  In  dit  rapport  is  uitgegaan  van  een  drukprofiel  dat  in  de  literatuur 
gegeven  wordt.  De  resultaten  van  het  model  werden  vergeleken  met  die  van  een 
empirische  formule  die  wordt  gebruikt  bij  de  damptesten  van  NBC-beschermende 
kledingmaterialen.  Bij  de  experimentele  condities  waarbij  de  empirische  formule 
werd  bepaald,  bleken  de  resultaten  een  goede  overeenkomst  te  vertonen.  Door 
middel  van  een  parameterstudie  werd  de  invloed  van  diverse  parameters  in  het 
model  onderzocht.  Indien  de  dikte  van  de  luchtlaag  onder  de  kleding  afneemt,  of  de 
diameter  van  de  cilinder  toeneemt,  kan  dit  leiden  tot  een  afhame  van  de  luchtsnel- 
heid  door  het  kledingmateriaal.  De  grootste  effecten  worden  gevonden  bij  mate- 
rialen  met  een  lage  luchtweerstand. 
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Romaans 

g 

Valversnelling 

(m/s2) 

h 

Hoogte 

(m) 

K 

Constante 

(-) 

Kb 

Berekende  constante 

(mmH20/(cm/s)/ (m/s» 

Ke 

Empirische  constante 

(mmH20/(cm/s)/ (m/s)) 

L 

Lengte  van  de  cilinder 

(m) 

Po 

Druk  ver  van  de  cilinder 

(Pa) 

Pin 

Druk  in  de  luchtlaag  onder  het  kledingmateriaal 

(Pa) 

Puit 

Druk  buiten  het  kledingmateriaal 

(Pa) 

R 

Straal  van  de  cilinder 

(m) 

Re 

Reynolds-getal 

(-) 

AR 

Halve  dikte  van  de  luchtlaag  tussen  wand  en  kledingmateriaal 

(m) 

Rinat 

Luchtweerstand  van  het  kledingmateriaal 

(mmH20/(cm/s)) 

VO 

Windsnelheid 

(m/s) 

vmat 

Luchtsnelheid  door  het  kledingmateriaal 

(m/s) 

X 

Plaats  langs  de  cilinder  ten  opzichte  van  het  stagnatiepunt 

(m) 

Grieks 

r 

Luchtpermeabiliteit 

(m/(Pa  s)) 

(|>V 

Debiet  in  de  luchtlaag 

(m3/s) 

T) 

Dynamische  viscositeit 

(kg/(m  s)) 

0 

Hoek  ten  opzichte  van  de  windrichting 

(°) 

Pi 

Dichtheid  lucht 

(kg/m3) 

Pw 

Dichtheid  water 

(kg/m3) 

X 

Schuifspanning 

(Pa) 
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1  Inleiding 


In  eerder  onderzoek  is  een  model  opgesteld,  dat  de  penetratie  van  damp  door 
NBC-beschermend  kledingmateriaal  beschrijft  [1-4].  Om  tot  een  volledige  model¬ 
iering  van  de  bescherming  van  een  totaal  kledingsysteem  te  komen,  moeten  de 
volgende  processen  beschreven  worden: 

1  luchtstromingsprofiel  om,  door  en  onder  de  kleding; 

2  damppenetratie  door  het  kledingmateriaal; 

3  damptransport  naar  de  huid. 

In  de  literatuur  wordt  een  fysisch  model  beschreven  voor  processen  die  de  penetra¬ 
tie  en  de  huiddepositie  van  aerosolen  door  permeabele  kledingsystemen  bepalen 
[5].  Het  lichaam(sdeel)  wordt  hierbij  gekarakteriseerd  door  een  cilinder  die  op  een 
bepaalde  afstand  door  kledingmateriaal  is  omspannen.  Met  behulp  van  dit  model 
kan  de  luchtstroming  door  het  materiaal,  de  penetratie  van  het  aerosol  door  het 
materiaal,  en  de  depositie  van  het  aerosol  op  het  cilinderoppervlak  worden  bere- 
kend. 

Voor  de  modeliering  van  de  dampbescherming  van  NBC-kleding  is  voor  een 
dergelijke  opzet  gekozen.  Allereerst  is  een  model  ontwikkeld  dat  het  luchtsto- 
mingsprofiel  om,  door  en  onder  de  kleding  beschrijft  (1).  In  combinatie  met  het 
model  voor  de  penetratie  van  damp  door  het  kledingmateriaal  (2)  en  het  nog  te 
ontwikkelen  model  voor  het  damptransport  naar  de  huid  (3),  kan  dit  leiden  tot  een 
model  voor  de  bescherming  die  door  een  compleet  kledingsysteem  geboden  wordt. 
In  dit  rapport  wordt  een  beschrijving  gegeven  van  het  model  voor  het  luchtstro¬ 
mingsprofiel  om,  door  en  onder  de  kleding.  Tevens  worden  de  resultaten  gegeven 
van  een  parameterstudie,  waarbij  de  invloed  van  de  diameter  van  de  cilinder,  de 
afstand  tussen  de  kleding  en  het  cilinderoppervlak,  en  de  luchtweerstand  van  het 
kledingmateriaal  op  het  stromingsprofiel  werd  bepaald. 
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2  Theorie 


2.1  Procesbeschrijving 

De  luchtsnelheid  door  het  kledingmateriaal  is  afhankelijk  van  de  aanstroomrich- 
ting  van  de  wind.  Aan  de  voorkant  van  het  lichaamsdeel  zal  de  lucht  naar  binnen 
stromen  en  aan  de  achterkant  zal  de  lucht  weer  naar  buiten  stromen.  In  principe  is 
op  elke  plaats  de  luchtsnelheid  door  het  kledingmateriaal  anders.  Het  is  mogelijk 
deze  stroming  te  modelleren  door  aan  te  nemen  dat  het  desbetreffende  lichaams¬ 
deel  een  cilinder  is  [5].  In  figuur  1  is  een  schematische  weergave  van  deze  cilinder 
met  het  stromingsprofiel  daaromheen  gegeven. 


Figuur  1:  Het  stromingsprofiel  van  lucht  om  een  cilinder,  bekleed  met  kledingmate¬ 

riaal. 

De  cilinder  is  bekleed  met  kledingmateriaal.  Dit  kledingmateriaal  zit  op  een  be- 
paalde  gedefinieerde  afstand  van  de  cilinder  (2AR),  zodat  er  een  luchtlaag  tussen 
het  kledingmateriaal  en  de  cilinder  zit.  Op  punt  C  raakt  de  luchtstroom  de  kleding 
loodrecht.  Dit  punt  wordt  het  stagnatiepunt  genoemd.  Op  andere  plaatsen  (A), 
onder  een  hoek  0  ten  opzichte  van  het  stagnatiepunt,  zal  de  luchtstroom  de  kleding 
niet  loodrecht  raken.  Als  gevolg  hiervan  zal  op  elke  plek  om  de  kleding  heen  een 
andere  buitendruk  heersen.  Deze  buitendruk  moet  bekend  zijn  om  het  totale  stro¬ 
mingsprofiel  te  berekenen. 

Ook  de  druk  onder  de  kleding  zal  afhankelijk  zijn  van  de  plaats.  Als  gevolg  van 
deze  plaatsafhankelijkheid  zal  lucht  in  de  luchtlaag  tussen  het  kledingmateriaal  en 
de  cilinder  door  stromen.  Deze  stroming  is  te  benaderen  als  luchtstroming  tussen 
twee  platen  door. 

Omdat  de  druk  onder  de  kleding  en  buiten  de  kleding  niet  gelijk  is,  zal  ook  lucht 
door  het  doek  gaan  stromen.  Deze  luchtstroom  is  te  berekenen  als  beide  drukken 
bekend  zijn. 
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2.2  Turbulentie  criteria 

Het  stromingsprofiel  om  een  cilinder  met  kledingmateriaal  is  sterk  afhankelijk  van 
de  mate  van  turbulentie.  Turbulentie  is  te  beschrijven  als  de  verhouding  tussen  de 
convectieve  kracht  (impuls)  en  de  viscositeit  (traagheid).  Een  bepaald  medium, 
bijvoorbeeld  lucht,  heeft  een  dichtheid  p.  Dit  medium  stroomt  met  een  bepaalde 
snelheid,  v,  tegen  een  bolletje  of  een  cilinder  met  een  straal,  R,  aan.  Als  gevolg  van 
deze  aanstroming  zullen  wervels  kunnen  ontstaan  in  de  lucht  achter  de  cilinder. 
Deze  wervels  zullen  echter  minder  optreden,  als  de  viscositeit  van  het  medium,  tj, 
hoger  is.  In  de  praktijk  is  gebleken,  dat  de  verhouding  tussen  deze  vier  variabelen 
een  goede  maat  vormt  voor  turbulentie: 


Re  = 


2pvR 


V 


(1) 


Deze  verhouding  wordt  het  kengetal  van  Reynolds  (Re)  genoemd.  Als  Re  laag  is, 
heet  de  stroming  laminair.  Wordt  Re  hoger  dan  2000  a  2400,  dan  zal  de  stroming 
turbulent  worden  en  zullen  er  wervelingen  gaan  plaatsvinden. 


2.3  Het  snelheidsprofiel  om,  door  en  onder  de  kleding 


2.3.1  De  luchtsnelheid  om  de  kleding 

Voor  het  berekenen  van  de  luchtsnelheid  door  de  kleding,  moet  het  drukprofiel 
buiten  de  cilinder,  bedekt  met  kleding,  bekend  zijn.  Dat  drukprofiel  is  afhankelijk 
van  de  mate  van  turbulentie  om  de  kleding. 

In  het  laminaire  geval,  zal  de  stroming  er  uitzien,  zoals  weergegeven  in  figuur  2. 


Figuur  2:  Stromingsprofiel  van  een  gasstroom  om  een  cilinder  bij  lage  Reynolds. 

Het  drukprofiel  om  een  cilinder  is  bekend  indien  de  luchtstroom  laminair  is.  In  dat 
geval  kan  dit  profiel  analytisch  berekend  worden  (bijlage  A).  Het  resultaat  is: 

Puit0  =  p0  +  ^Plv0  t  ~4sin2  G) 


(2) 
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Waarbij  Puit,  0  de  plaatselijke  druk  is  op  een  plaats,  die  een  hoek  0  verwijderd  is 
van  het  stagnatiepunt  (punt  C  in  figuur  2).  Deze  druk  is  afhankelijk  van  de  buiten- 
druk  voor  de  cilinder,  Pq,  de  windsnelheid,  vq,  en  de  dichtheid  van  de  lucht,  p\. 

Indien  de  stroming  om  de  cilinder  turbulent  is,  zullen  wervels  achter  de  cilinder 
ontstaan.  De  druk  om  de  cilinder  is  dan  niet  te  berekenen  en  dient  experimented  te 
worden  bepaald. 


2.3.2  De  luchtsnelheid  door  de  kleding 

Als  gevolg  van  het  drukverschil,  AP,  tussen  de  plaatselijke  druk  buiten  het  kle- 
dingmateriaal  en  de  druk  in  de  luchtlaag  tussen  het  kledingmateriaal  en  de  cilinder, 
zal  de  lucht  door  het  kledingmateriaal  heen  kunnen  stromen.  In  bijlage  B  is  afge- 
leid,  dat  de  snelheid  waarmee  de  lucht  door  het  kledingmateriaal  stroomt,  vmat, 
berekend  kan  worden  met: 


10AP 

R matPwS 


(3) 


Waarbij  Rmat  de  luchtweerstand  van  het  kledingmateriaal  is,  pw  de  dichtheid  van 
water  en  g  de  valversnelling.  Om  deze  snelheid  te  berekenen  moeten  dus  zowel  de 
druk  buiten  de  kleding  als  de  druk  onder  de  kleding  bekend  zijn.  Uit  de  plaatselijke 
binnen-  en  buitendruk  kan  dan  de  plaatselijke  luchtsnelheid  door  het  kledingmate¬ 
riaal  berekend  worden  met  behulp  van  vergelijking  (3). 


2.3.3  De  luchtsnelheid  onder  de  kleding 

In  bijlage  C  is  afgeleid,  dat  het  debiet  van  de  stroming  tussen  de  kleding  en  de 
cilinder,  <|>v,  als  volgt  berekend  kan  worden,  mits  de  stroming  onder  de  kleding 
laminair  is: 


3  7 


dpA  y^y 

dd  g  R 


(4) 


De  dikte  van  de  luchtlaag  onder  de  kleding  is  hierbij  2AR,  de  lengte  van  de  cilinder 
is  L  en  de  straal  is  R.  De  viscositeit  van  de  lucht  is  gelijk  aan  tv  Het  debiet  veran- 
dert  als  functie  van  de  hoek,  0.  De  verandering  van  het  debiet  is  gerelateerd  aan  de 
luchtsnelheid  door  het  kledingmateriaal,  vmat.  In  formulevorm  kan  dat  worden 
weergegeven  als  [5]: 

d<pv  =  v^Ldx  =  v^LRdG  (5) 


Waarbij  x  de  plaats  op  de  cilinder  is  ten  opzichte  van  het  stagnatiepunt.  Als  verge¬ 
lijking  (3)  hierin  ingevuld  wordt,  wordt  gevonden  dat: 


d<pv(9)= 


10LRAP 

R matPwS 


de 


(6) 


Vervolgens  kan  vergelijking  (4)  hierin  worden  ingevuld.  Dit  geeft: 
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2L 

r 

d2p ■ 

u  rin 

\ 

(ptf 

10LRAP 

10LR 

3jj 

\ 

dd2 

R 

Rmat  PwS 

R mat  Pwg 

(7) 


Deze  vergelijking  beschrijft  het  drukprofiel  tussen  de  kleding  en  de  cilinder  als 
functie  van  de  plaats.  Met  behulp  van  deze  vergelijking  kan  de  plaatselijke  druk 
tussen  het  kledingmateriaal  en  de  cilinder  berekend  worden.  Gecombineerd  met 
vergelijking  (4)  kan  het  snelheidsprofiel  onder  de  kleding  berekend  worden. 
Overigens  dient  opgemerkt  te  worden,  dat  deze  vergelijking  niet  meer  opgaat  bij 
het  zogenaamde  ‘close  contact’,  waarbij  de  afstand  tussen  de  kleding  en  de  huid 
nul  is. 


2.3.3.1  Laminaire  luchtsnelheid 

Als  de  luchtsnelheid  buiten  de  kleding  laminair  is,  is  het  mogelijk  vergelijking  (7) 
analytisch  op  te  lossen  (bijlage  D).  De  luchtsnelheid  door  het  kledingmateriaal  is  in 
dat  geval: 


v  = 


Pi 

PwS 


1  + 


lo(l-2sin2  o) 

30  tj  R2 
8  pwg  (AR  y  Rmat 


vo 


R 


mat 


(8) 


of  lokaal,  bij  een  hoek  0  ten  opzichte  van  het  stagnatiepunt: 


K mat 


Kh 


v2o 


R 


mat 


(9) 


Waarin  Kb  een  constante  is,  die  afhangt  van  de  straal  van  de  cilinder,  de  dikte  van 
de  luchtlaag,  de  luchtweerstand  van  het  kledingmateriaal  en  de  hoek  ten  opzichte 
van  het  stagnatiepunt. 


2.3.3.2  Turbulente  luchtsnelheid 

Indien  de  stroming  om  de  kleding  turbulent  is,  kan  het  drukprofiel  niet  berekend 
worden.  Dit  profiel  zal  dan  moeten  volgen  uit  metingen.  In  dat  geval  is  vergelijking 
(7)  niet  analytisch  oplosbaar.  Het  is  wel  mogelijk  om  met  behulp  van  numerieke 
berekeningen  dit  profiel  te  berekenen.  Hiervoor  is  krachtige  software  nodig.  Aan- 
gezien  deze  software  niet  voorhanden  is,  is  gekozen  om  gebruik  te  maken  van 
gemeten  drukprofielen. 

Het  is  wel  mogelijk  om  een  schatting  te  maken  van  de  gemiddelde  snelheid  door 
het  kledingmateriaal  (bijlage  E).  In  dat  geval  kan  voor  de  gemiddelde  luchtsnelheid 
door  het  kledingmateriaal  gevonden  worden: 


v  _  5  Pi  vq 

PwS  Rmat 


(10) 


of  weer: 
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vmat  ~  Kb 


v0 

R 


mat 


(9) 


Waarbij  alleen  de  waarde  van  Kb  verschilt  ten  opzichte  van  het  laminaire  geval. 
Deze  snelheid  is  het  gemiddelde  over  alle  plaatsen  waar  de  lucht  naar  binnen 
stroomt. 


2.4  Empirische  relatie  voor  het  berekenen  van  de  luchtsnelheid 


Bij  het  bepalen  van  de  in  te  stellen  luchtsnelheid  bij  damptesten  van  NBC- 
beschermende  kleding,  wordt  een  empirische  formule  gebruikt  om  de  luchtsnelheid 
door  het  kledingmateriaal  te  voorspellen  als  functie  van  de  windsnelheid  en  de 
luchtweerstand  van  het  kledingmateriaal  [6-8]: 


vmat 


(11) 


Waarbij  Ke  een  empirische  constante  is:  Ke  =  5,7  *  10‘4  (mmH20/(cm/s))/(m/s). 
Deze  empirische  formule  is  gevonden  bij  de  volgende  omstandigheden  [8]: 

R  =  0,025  m; 

2AR  =  2  *  lO'3 m; 

Rmat  =  2,66  mmH20/(cm/s); 
vo  =  9,1  m/s. 

Deze  luchtsnelheid  is  een  gemiddelde  over  het  gedeelte  van  het  kledingmateriaal, 
waar  de  luchtstroom  naar  binnen  gaat. 


Vergelijking  (9)  toont  een  grote  gelijkenis  met  de  empirische  vergelijking  (11). 
Zowel  vergelijking  (8)  als  (10)  geven  een  analytische  uitdrukking  voor  de  empiri¬ 
sche  constante  K*. 


Deze  empirische  formule  is  vergeleken  met  het  model. 
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3  Berekeningen 


3.1  Experimentele  omstandigheden 

Om  drukprofielen  te  berekenen,  moeten  verschillende  constanten  en  variabelen 
bekend  zijn.  Er  wordt  uitgegaan  van  standaardomstandigheden  (gedefinieerde 
constante  druk  en  temperatuur  etcetera  [9]),  zodat  geldt: 
r\  =  18,1  *10-6  Pas; 

pi  =1,182  kg/m3; 
pw  =  1000  kg/m3; 
g  =9,81  m/s2. 

Voor  de  straal  van  de  cilinder  is  gekozen  voor  R  =  0,05  m,  wat  overeenkomt  met 
de  straal  van  een  arm.  De  dikte  van  de  luchtlaag  tussen  de  kleding  en  de  huid  is 
gekozen  als  2AR  =  2*  10'3  m. 

De  luchtweerstanden  van  een  drietal  karakteristieke  typen  NBC-kleding  zijn: 
USM  :  Rmat  =  2,21  mmH20/(cm/s); 

M82  :  Rmat  =  0,91  mmH20/(cm/s); 

OOA  :  Rmat  =0,16  mmH20/(cm/s). 

Voor  de  windsnelheid  werd  gekozen  voor  1  m/s  en  5  m/s.  Deze  waarden  zijn 
gebruikelijk  voor  NBC-kledingmateriaaltesten  [7]. 


3.2  Luchtsnelheid 


3.2.1  Laminaire  stroming 

Indien  de  stroming  om  de  kleding  laminair  is,  beschrijft  formule  (8)  de  relatie 
tussen  de  druk  en  de  snelheid  als  functie  van  de  aanstroomhoek  (Re  <  50).  Bij  een 
snelheid  van  vo  =  1  m/s  is  de  stroming  om  de  kleding  duidelijk  turbulent 
(Re  =  6530).  Om  een  indicatie  te  krijgen,  werd  toch  op  basis  van  laminaire  verge- 
lijkingen  de  drukverdeling  buiten  en  onder  de  kleding  berekend  (figuur  3). 
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Figuur  3:  De  relatieve  drukverdeling  buiten  en  onder  de  kleding  bij  een  windsnelheid 

van  1  m/s  (Rmat  ~  0,16  mmfyOffcmfs)). 

De  weergegeven  drukken  zijn  relatieve  drukken  ten  opzichte  van  de  druk  ver  van 
de  cilinder.  Om  die  reden  kunnen  ook  negatieve  drukken  gevonden  worden.  Uit  dit 
drukprofiel  kan  de  plaatselijke  snelheid  door  het  kledingmateriaal  berekend  worden 
(figuur  4). 


Figuur  4:  De  snelheidsverdeling  door  de  kleding  bij  een  windsnelheid  van  1  m/s, 

waarbij  een  positieve  waarde  stoat  voor  een  luchtsnelheid  naar  binnen  en 
een  negatieve  waarde  voor  een  luchtsnelheid  naar  buiten  (Rmat  -  0,16  mm 
H20f(cm/s)). 

Voor  plaatsen  waar  de  buitendruk  hoger  is  dan  de  binnendruk,  zal  de  lucht  naar 
binnen  stromen,  terwijl  in  het  tegengestelde  geval  de  lucht  naar  buiten  stroomt. 
Vergeleken  met  de  buitendruk,  is  de  binnendruk  bijna  constant.  Omdat  de  snelheid 
door  het  kledingmateriaal  schaalt  met  het  drukverschil  tussen  de  binnen-  en  de 
buitendruk,  schaalt  in  dit  geval  de  snelheid  met  de  buitendruk. 
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3.2.2  Turbulente  stroming 

Bij  turbulente  stroming  om  de  kleding  is  het  drukprofiel  om  de  cilinder  anders.  De 
aanpak  van  turbulente  stroming  gaat  uit  van  experimented  gemeten  buitendrukken. 
In  de  literatuur  worden  gemeten  buitendrukken  gegeven  als  functie  van  de  aan- 
stroomhoek  om  een  cilinder,  waarbij  Reynolds  gelijk  is  aan  1,9  *  105  [10].  Dit  is 
bijvoorbeeld  het  geval  bij  een  cilinder  met  een  straal  van  0,1  m  bij  een  windsnel- 
heid  van  14,5  m/s  (de  overige  parameters  zijn  gelijk  gesteld  aan  de  vorige  bereke- 
ningen).  Volgens  [11]  blijft  dit  drukprofiel  qua  vorm  bijna  constant  als  de  mate  van 
turbulentie  varieert.  De  diameter  van  de  cilinder  is  dus  voor  dit  profiel  vrij  onbe- 
langrijk.  Om  het  drukprofiel  te  krijgen  bij  een  andere  windsnelheid,  moet  dan 
alleen  gecorrigeerd  worden  ten  opzichte  van  het  stagnatiepunt  (daar  is  de  druk 
gelijk  aan: 

1  2 

~Plvmat’ 

zodat  het  totale  profiel  gecorrigeerd  kan  worden).  Volgens  [1 1]  is  deze  werkwijze 
gerechtvaardigd,  aangezien  de  metingen  aantoonden  dat  het  drukprofiel  om  een 
cilinder  hetzelfde  patroon  vertoont  tussen  5  m/s  en  17  m/s. 

In  figuur  5  zijn  de  gemeten  drukken  weergegeven,  samen  met  de  met  vergelijking 
(7)  berekende  binnendruk. 


Figuur  5:  De  gemeten  relatieve  buitendruk  [10]  en  de  daaruit  berekende  relatieve 

binnendruk  als  functie  van  de  aanstroomhoek  bij  een  windsnelheid  van 
14,5  mis  en  een  cilinderdiameter  van  0,1  m. 

Uit  deze  twee  drukken  kan  vervolgens  de  luchtsnelheid  door  het  kledingmateriaal 
berekend  worden  met  behulp  van  vergelijking  (3).  Deze  luchtsnelheid  is  weergege¬ 
ven  in  figuur  6.  Dit  is  gedaan  voor  drie  verschillende  materialen  met  verschillende 
luchtweerstand. 
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Figuur  6:  De  luchtsnelheid  door  kledingmateriaal  als  Junctie  van  de  aanstroomhoek , 

volgend  uit  drukmetingen  [10]  bij  varierende  luchtweerstand  (vq  =  14 >5  mis 
enR  ~  0,1  mj 

Als  dit  profiel  gebruikt  wordt  voor  een  windsnelheid  van  5  m/s,  wordt  een  analoog 
profiel  gevonden  (figuur  7). 


Figuur  7:  De  luchtsnelheid  door  kledingmateriaal  als  functie  van  de  aanstroomhoek, 

volgend  uit  drukmetingen  [10]  bij  varierende  luchtweerstand  (vq  =  5  mis  en 
R  =  0,1  m). 
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De  maximale  snelheid  van  de  lucht  in  de  luchtlaag  tussen  het  kledingmateriaal  en 
de  cilinder  werd  als  functie  van  de  aanstroomhoek  bepaald  volgens  bijlage  B.  De 
resultaten  zijn  te  zien  in  figuren  8  en  9. 


Figuur  8:  De  luchtsnelheid  tussen  kledingmateriaal  en  de  cilinder  als  functie  van  de 

aanstroomhoek :,  bij  varierende  luchtweerstand  (v  =  14,5  m/s  enR  =  0,1  m). 
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Figuur  9:  De  luchtsnelheid  tussen  kledingmateriaal  en  de  cilinder  als  functie  van  de 

aanstroomhoek,  bij  varierende  luchtweerstand  (v  =  5  m/s  enR  =  0,1  m). 

Opvallend  is,  dat  bij  hoge  windsnelheid  en  lage  luchtweerstand,  de  luchtsnelheid 
bij  een  bepaalde  hoek  bijna  even  groot  is  als  de  windsnelheid.  Het  is  echter  de 
vraag  of  de  gebruikte  vergelijkingen  in  dat  regime  nog  geldig  zijn.  Zo  is  de  stro- 
ming  tussen  de  cilinder  en  het  kledingmateriaal  nog  wel  laminair,  maar  al  vrij  groot 
(Re  =  1894).  Indien  die  stroming  turbulent  wordt,  geldt  vergelijking  (7)  niet  meer. 
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Ook  de  luchtstroora  door  het  kledingmateriaal  is  vrij  groot.  Dit  kan  tot  gevolg 
hebben,  dat  de  luchtsnelheid  door  het  kledingmateriaal  niet  meer  lineair  met  het 
drukverschil  schaalt.  In  bepaalde  gevallen  kan  een  kwadratisch  verband  gevonden 
worden  [12, 13]. 


3.3  Parameterstudie 

Voor  turbulente  stroming  is  de  invloed  van  verschillende  parameters  op  het  stro- 
mingsprofiel  door  de  kleding  bestudeerd.  Het  stromingsprofiel  wordt  hierbij  geka- 
rakteriseerd  door  Kb  (paragraaf  2.3.3.2)  en  de  maximale  hoek  waarbij  de  lucht- 
stroom  door  het  kledingmateriaal  inwaarts  gericht  is  (bijlage  D).  Kb  wordt  bere- 
kend  uit  de  gemiddelde  luchtsnelheid  door  het  kledingmateriaal  op  de  plaatsen 
waar  deze  naar  binnen  toe  gericht  is.  De  empirische  waarde  voor  deze  parameter  is 
gelijk  aan  Ke  =  5,7  *  10‘4  (mmH20/(cm/s))/(m/s)  (paragraaf  2.4)  [8]. 

De  straal  van  de  cilinder  werd  gevarieerd  van  0,025  m  (een  vinger)  tot  0,3  m  (de 
romp).  De  dikte  van  de  luchtlaag  onder  de  kleding  werd  gevarieerd  van  0,001  m  tot 
0,01  m.  Dit  is  gedaan  voor  de  drie  kledingmaterialen  met  verschillende  luchtweer- 
stand,  zoals  genoemd  in  paragraaf  3.1. 

3.3.1  De  dikte  van  de  luchtlaag  onder  de  kleding 

De  invloed  van  dikte  van  de  luchtlaag  op  Kb  en  0  is  te  zien  in  figuren  10  en  1 1. 
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Figuur  10:  De  invloed  van  de  dikte  van  de  luchtlaag  op  de  constante  Kb  bij  verschillende 
typen  kledingmateriaal  (R  =  0,1  m). 
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Figuur  11:  De  invloed  van  de  dikte  van  de  luchtlaag  op  de  maximale  hoek  van  inwaartse 
luchtstroom  bij  verschillende  typen  kledingmateriaal  (R  =  0,1  m). 

Als  de  luchtlaag  dikker  is  dan  5  mm,  zijn  zowel  Kb  als  0  vrijwel  constant.  Kb  is 
circa  6,2  *  10'4  (mmH20/(cm/s)/(m/s)).  Deze  waarde  blijkt  goed  overeen  te  komen 
met  de  experimented  bepaalde  waarde  van  Ke  =  5,7  *  10'4  (mmH20/(cm/s)/(m/s)) 
(paragraaf  2.4).  Als  de  luchtlaag  dunner  wordt,  kan  zowel  Kb,  als  de  hoek  kleiner 
worden.  Wanneer  dit  gebeurt,  is  afhankelijk  van  de  luchtweerstand  van  het  kle¬ 
dingmateriaal.  Het  effect  is  groter  bij  materialen  met  een  lage  luchtweerstand.  Het 
is  opvallend,  dat  de  maximale  instroomhoek  afneemt  bij  afnemende  dikte  van  de 
luchtlaag,  maar  bij  verdere  verlaging  weer  gaat  toenemen.  Een  mogelijke  oorzaak 
is,  dat  de  snelheid  door  het  kledingmateriaal  bij  lage  luchtweerstanden  niet  meer 
lineair  verloopt  met  het  drukverschil  [12].  Het  is  echter  ook  mogelijk,  dat  de  ge- 
bruikte  vergelijkingen  voor  het  berekenen  van  het  drukprofiel  onder  de  kleding  bij 
die  situatie  niet  meer  toepasbaar  zijn.  Om  een  uitspraak  hierover  te  doen,  moeten 
drukmetingen  verricht  worden  aan  dit  systeem.  Dit  kan  gebeuren  door  zowel  de 
binnen-  als  de  buitendruk  te  meten  bij  een  met  lucht  aangestroomde  cilinder,  be- 
kleed  met  kledingmateriaal.  Met  behulp  van  de  resultaten  van  die  metingen  kunnen 
de  gebruikte  vergelijkingen  geverifieerd  worden. 

Als  de  dikte  van  de  luchtlaag  gelijk  wordt  aan  0  (‘close  contact’)  zullen  de  ge¬ 
bruikte  vergelijkingen  in  ieder  geval  niet  meer  opgaan.  Daar  zal  een  ander  model 
voor  ontwikkeld  moeten  worden. 

3.3.2  De  straal  van  de  cilinder 

In  figuren  12  en  13  is  het  effect  van  de  straal  van  de  cilinder  op  Kb  en  0  te  zien. 
Hierbij  is  aangenomen,  dat  het  drukprofiel  om  de  cilinder  niet  verandert,  hoewel 
Reynolds  verandert  met  de  straal.  In  werkelijkheid  zal  dit  profiel  dus  iets  kunnen 
afwijken.  In  welke  mate  zal  uit  drukmetingen  moeten  volgen. 
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Figuur  12:  De  invloed  van  de  straal  van  de  cilinder  op  Kb  bij  verschillende  typen  kle- 
dingmateriaal  (2AR  -2  *  10~^  m). 


R(m) 


Figuur  13:  De  invloed  van  de  straal  van  de  cilinder  op  de  maximale  hoek  van  inwaartse 
luchtstroom  bij  verschillende  typen  kledingmateriaal  (2AR  =  2  *  10'^  m). 

Een  grotere  cilinderstraal  leidt  tot  verkleining  van  Kb*  Verder  heeft  ook  hier  verla- 
ging  van  de  luchtweerstand  van  kledingmateriaal  als  effect,  dat  de  verschillende 
gevolgen  eerder  merkbaar  worden.  0  Vertoont  een  iets  ander  gedrag.  Bij  hoge 
luchtweerstand  zal  vergroting  van  de  straal  leiden  tot  afname  van  0.  Bij  lagere 
luchtweerstand  is  dit  in  eerste  instantie  ook  het  geval,  maar  bij  een  bepaalde  straal 
gaat  deze  hoek  weer  toenemen.  Dit  komt  doordat  de  luchtweerstand  door  het  doek 
kleiner  wordt  ten  opzichte  van  de  weerstand  van  de  stroming  om  het  doek  (de 
cilinder). 


3.3.3  De  luchtweerstand  van  de  kleding 

De  invloed  van  de  luchtweerstand  van  het  kledingmateriaal  op  Kb  en  0  wordt  in 
figuur  14  gegeven. 
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Figuur  14:  De  invloed  van  de  luchtweerstand  van  kledingmateriaal  op  Kb  en  op  de 

maximale  hoek  van  inwaartse  luchtstroom  (v^t  =  $  m^s>  2AR  =  2  *  10'^  m, 
R  =  OJ  m). 

Verlaging  van  de  luchtweerstand  geeft  in  het  algemeen  een  verlaging  van  Kb,  maar 
als  de  luchtweerstand  groter  is  dan  circa  1  mmH20/(cm/s)  is  dit  effect  klein.  Ook 
de  maximaal  inwaartse  hoek  neemt  tot  een  bepaald  punt  af.  Bij  zeer  lage  lucht- 
weerstanden  (Rmat  <0,16  mmH20/(cm/s))  neemt  deze  hoek  weer  toe,  maar  dit 
effect  zou  een  numeriek  probleem  kunnen  zijn,  omdat  vergelijking  (7)  dan  mis- 
schien  niet  meer  toepasbaar  is.  Drukmetingen  kunnen  hierover  uitsluitsel  geven. 
Figuur  14  verandert  niet  als  functie  van  de  windsnelheid,  omdat  het  drukprofiel  om 
de  cilinder  niet  verandert  bij  verschillende  windsnelheden. 


3.3.4  Vergelijking  berekende  met  gemeten  waarde 
Uit  het  voorgaande  is  gebleken  dat  de  mate  van  overeenkomst  tussen  Kb  en  Ke 
afhankelijk  is  van  de  straal  van  de  cilinder,  R,  de  dikte  van  de  luchtlaag,  2AR,  en  de 
luchtweerstand  van  het  materiaal,  Rmat*  De  procentuele  afwijking  (verlaging)  van 
Kb  ten  opzichte  van  de  experimentele  waarde,  Kg,  als  functie  van  R  en  2AR  is 
weergegeven  in  figuren  15,  16  en  17  voor  drie  verschillende  typen  kledingmateri- 
alen,  waarbij  het  percentage  volgt  uit: 


Percentage  =  100* 


(12) 


Ook  de  omstandigheden,  waarbij  0%  afwijking  ten  opzichte  van  de  experimentele 
waarde  gevonden  wordt,  staan  in  deze  figuren  gegeven.  De  experimentele  omstan¬ 
digheden,  waarbij  Ke  bepaald  is,  zijn:  2AR  =  2*10  3  m  en  R  =  0,025  m  (para- 
graaf  2.4).  De  waarden  van  Kb  zijn  berekend  met  behulp  van  de  fit,  zoals  beschre- 
ven  in  bijlage  F.  Deze  fit  geeft  een  correlatiecoefficient  van  0,991  met  de  uit  het 
model  berekende  waarden. 


2AR  (m) 
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Figuur  17:  Procentuele  verlaging  van  Kb  ten  opzichte  van  de  experimentele  waarde,  Ke, 
als  functie  van  2AR  en  R,  bij  kledingmateriaal  OOA. 

Uit  de  figuren  is  te  zien,  dat  vergroting  van  de  cilinderstraal,  R,  of  verkleining  van 
de  dikte  van  de  luchtlaag,  2AR,  een  verlaging  van  de  waarde  van  Kb  tot  gevolg 
heeft.  Verder  is  te  zien,  dat  deze  effecten  groter  zijn  bij  kledingmateriaal  wat  een 
lagere  luchtweerstand  heeft. 

Aan  de  andere  kant  zal  bij  toename  van  2AR  en  afname  van  R  de  berekende  waar¬ 
de,  Kb,  naderen  tot  een  limietwaarde  5,81  *  10'4  (mmH20/(cm/s)/(m/s)).  Dit  is 
1,9%  boven  de  gemeten  waarde,  Kg.  Indien  de  limietwaarde  niet  uit  de  bovenge- 
noemde  fit  berekend  wordt,  maar  uit  het  model  moet  volgen,  zal  deze  niet  berekend 
kunnen  worden,  maar  zal  benaderd  moeten  worden.  In  paragraaf  3.3.1  werd  op  die 
manier  een  waarde  van  circa  6,2  *  10"4  (mmH20/(cm/s)/(m/s))  gevonden,  hetgeen 
8,8%  boven  de  experimentele  waarde  ligt. 

De  waarde  van  Ke  komt  dus  goed  overeen  met  de  maximale  waarde  van  Kb-  Aan- 
gezien  de  damptest  voor  NBC-beschermende  kledingmaterialen  wordt  uitgevoerd 
bij  een  luchtsnelheid  door  het  kledingmateriaal  die  evenredig  is  met  Kg,  kan  deze 
testopzet  als  ‘worst  case’-scenario’  worden  beschouwd. 
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4  Conclusies 


Er  is  een  model  ontwikkeld,  dat  de  luchtstroming  om,  door  en  onder  de  kleding  om 
een  lichaamsdeel  kan  berekenen.  Het  lichaamsdeel  wordt  hierbij  gekarakteriseerd 
door  een  cilinder,  die  op  een  bepaalde  afstand  door  kledingmateriaal  is  omspannen. 

Als  de  stroming  om  de  kleding  laminair  is,  zijn  stromingsprofielen  om,  door  en 
onder  de  kleding  analytisch  te  berekenen.  In  de  praktijk  echter,  is  de  luchtstroom 
om  de  kleding  vaak  turbulent.  In  dat  geval  moet  het  drukprofiel  om  de  kleding 
bekend  zijn  om  het  stromingsprofiel  te  kunnen  berekenen.  Dit  drukprofiel  kan 
volgen  uit  numerieke  berekeningen  of  kan  experimenteel  bepaald  worden.  In  dit 
rapport  zijn  experimentele  gegevens  uit  de  literatuur  gebruikt. 

Bij  de  uitvoering  van  de  damptest  van  NBC-beschermende  kledingmaterialen 
wordt  gebruikgemaakt  van  een  empirische  formule  voor  de  berekening  van  de 
luchtsnelheid  door  het  materiaal.  Deze  formule  werd  vergeleken  met  de  uitkomsten 
van  berekeningen,  verkregen  met  behulp  van  het  model.  De  luchtsnelheid  door  de 
kleding  wordt  binnen  een  marge  van  10%  van  de  experimenteel  bepaalde  waarde 
voorspeld.  Dit  geldt  echter  alleen  bij  dezelfde  experimentele  omstandigheden,  als 
waarbij  de  empirische  formule  bepaald  is. 

Een  parameterstudie  toont  aan,  dat  de  luchtsnelheid  door  het  kledingmateriaal 
afneemt  bij  afnemende  dikte  van  de  luchtlaag  onder  de  kleding.  Het  grootste  effect 
wordt  gevonden  bij  kledingmateriaal  met  een  lage  luchtweerstand.  Vergroting  van 
de  diameter  van  de  cilinder  leidt  eveneens  tot  afname  van  de  luchtsnelheid  door  het 
kledingmateriaal.  Ook  hier  wordt  het  grootste  effect  gevonden  bij  kledingmateriaal 
met  lage  luchtweerstand. 

Ter  verbetering  van  het  model  zullen  drukprofielen  om  cilinders  bekleed  met 
kledingmateriaal  gemeten  worden  en/of  numeriek  berekend  worden  met  behulp 
van  een  ‘Computational  Fluid  Dynamics’-programma. 

Verder  zal  het  model  worden  uitgebreid  met  de  mogelijkheid  van  ‘close  contact’, 
waarbij  de  afstand  tussen  de  kleding  en  de  huid  gelijk  is  aan  nul. 

Het  model  zal  in  combinatie  met  het  reeds  bestaande  model  voor  de  damppenetra- 
tie  door  het  materiaal  als  uitgangspunt  dienen  voor  het  modelleren  van  het  damp- 
transport  naar  het  cilinderoppervlak. 
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Bijlage  A  Het  drukprofiel  om  een  cilinder  bij  laminaire 
stroming 


Als  een  gasstroom  tegen  een  cilinder  stroomt,  kan  een  stromingsprofiel  ontstaan, 
zoals  weergegeven  in  Figuur  A.l.  Voorwaarde  voor  het  ontstaan  van  een  dergelijk 
stromingsprofiel  is  dat  de  gassnelheid  erg  laag  is  (Reynolds  <  50).  In  dat  geval  zal 
het  stromingspatroon  achter  de  cilinder  precies  het  spiegelbeeld  zijn  van  voor  de 
cilinder.  Als  dat  niet  het  geval  is,  zal  de  gasstroom  aan  de  achterkant  van  de  cilin¬ 
der  eerder  loslaten,  zodat  het  stromingsprofiel  verandert.  Er  zullen  kleine  wervelin- 
gen  achter  de  cilinder  kunnen  ontstaan.  Het  grootste  deel  van  het  stromingsprofiel 
blijft  echter  over  het  totale  laminaire  gebied  gelijk. 


Figuur  A.l:  Stromingsprofiel  van  een  gasstroom  om  een  cilinder  bij  lage  Reynolds 
(Re  <  50). 

Op  het  punt  C  (het  stagnatiepunt)  is  de  gassnelheid  0.  Op  het  punt  A  stroomt  het 
gas  met  een  bepaalde  snelheid,  va,  langs.  Deze  snelheid  zal  in  dit  laminaire  geval 
gelijk  zijn  aan  [14]: 

vA=2v0sin6  (A.l) 

Langs  een  stroomlijn,  zoals  weergegeven  in  figuur  A.l,  kan  de  druk  op  elk  punt  op 
de  stroomlijn  berekend  worden  met  behulp  van  de  vergelijking  van  Bernoulli 
[A.1]: 

P  -  1  o  . . 

— I -  gh  +—v  =  const.  (A.2) 

P  2 

In  dit  geval  kan  aangenomen  worden,  dat  h=0.  Vergelijken  we  nu  de  situatie  op 
plek  0  en  op  plek  A,  dan  kan  voor  vergelijking  (A.2)  geschreven  worden: 

Pa  =Pq  +  j/A'o  ~jPva  (A.3) 
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Als  vergelijking  (A.  1 )  hierin  wordt  ingevuld,  levert  dit: 

PA=P0+^Pv0^~4  sin2  p)  (A.4) 

Met  behulp  van  deze  vergelijking  kan  dus  het  drukprofiel  om  een  cilinder  heen 
berekend  worden  in  het  geval  van  laminaire  stroming. 

[1 4]  Streeter,  V.L.  en  Wylie,  E.B., 

‘Fluid  mechanics;  first  SI  metric  edition’, 

McGraw-Hill,  Auckland,  2, 1987. 
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Bijlage  B  Luchtsnelheid  door  kledingmateriaal 


De  snelheid,  v  (in  m/s),  waarmee  de  lucht  door  het  kledingmateriaal  gaat,  is  lineair 
evenredig  met  het  drukverschil,  AP,  over  het  kledingmateriaal.  In  principe  kan 
kledingmateriaal  gezien  worden  als  een  filter.  De  wet  van  Darcy  beschrijft  de 
snelheid  door  een  filter  als  volgt: 

v  -K—AP  (B.l) 

V 

Waarin  K  een  constante  is,  1  de  dikte  van  het  filter  is  en  ri  de  viscositeit  van  het 
medium  is.  K,  1  en  T|  worden  in  het  geval  van  lucht  soms  gecombineerd  tot  een 
nieuwe  term,  T,  welke  de  luchtpermeabiliteit  (in  m/(Pa  s))  genoemd  wordt,  zodat: 

v  =  rAP  (B.2) 


De  reciproke  van  de  luchtpermeabiliteit  is  de  luchtweerstand.  Het  is  de  gewoonte 
deze  luchtweerstand  uit  te  drukken  in  mmH20/(cm/s).  Dit  levert  voor  de  lucht¬ 
weerstand,  Rinat,  op: 


R 


mat 


w 

^lPwg 


(B.3) 


Waarbij  pw  de  dichtheid  van  water  is  en  g  de  valversnelling  is.  Voor  de  snelheid 
door  het  kledingmateriaal  kan  dus  geschreven  worden: 


10AP 
R mat  Pw8 


(B.4) 


De  snelheid  van  lucht  door  kledingmateriaal  hangt  dus  af  van  de  luchtweerstand 
van  dat  kledingmateriaal. 
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Bijlage  C  Het  snelheidsprofiel  onder  kledingmateriaal 


In  eerste  instantie  kan  het  snelheidsprofiel  tussen  kledingmateriaal  en  een  cilinder 
benaderd  worden  met  het  snelheidsprofiel  tussen  twee  platen.  In  figuur  C.l  zijn 
twee  platen  weergegeven,  waar  gas  met  een  bepaald  debiet,  <|>v,  tussendoor  stroomt. 
De  afstand  tussen  de  twee  platen  is  gelijk  aan  2AR.  De  platen  steken  een  lengte  L 
uit  de  figuur  in  de  diepte. 


<t>v 


2AR 


Xi 


Figuur  C.l:  Luchtstroom  tussen  twee  platen  met  hierin  een  volume-elementje. 

Beschouwen  we  een  volume-elementje  met  een  dikte  2r,  een  lengte  fe-xj)  en  een 
diepte  L,  dan  kunnen  we  daar  een  krachtenbalans  over  berekenen.  De  krachten 
waar  rekening  mee  gehouden  moet  worden,  zijn  de  druk  die  voor  en  achter  het 
elementje  op  het  oppervlak  2rL  drukt,  en  de  schuifspanning,  x,  die  op  de  twee 
zijkanten  van  het  elementje  werkt.  In  formulevorm  ziet  deze  balans  er  uit  als: 

2rLPj  -  2rLP2  -2L(x2  -  jc,  >  =0  (C.  1) 


Als  de  afstand  tussen  xi  en  X2  klein  wordt  gekozen,  zodat  (X2-X1)  =  dx,  kan  deze 
vergelijking  geschreven  worden  als: 


De  schuifspanning  is  gerelateerd  aan  de  viscositeit,  rj,  van  het  medium  (de  lucht) 
en  het  snelheidsverschil  tussen  de  twee  lagen.  Voor  een  Newtons-medium  (zoals 
bijvoorbeeld  lucht  of  water)  geldt: 


(C.3) 


zodat: 


dv  _  r  r  dP  \ 
dr  Tjy  dx ) 


(C.4) 


Dus  voor  de  snelheid  op  punt  r,  tussen  de  platen,  kan  na  integratie  gevonden  wor¬ 
den: 


vr 


2rjy  dx  J 


(C.5) 
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Waarbij  A  een  integratieconstante  is.  Aan  de  rand  van  de  spleet,  op  r=AR,  is  de 
snelheid  vr  gelijk  aan  0.  Dit  levert  op,  dat: 


Vr 


(C.6) 


Dit  is  het  bekende  parabolische  snelheidsprofiel  van  laminaire  flow  (Poisseulle- 
stroming).  Hieruit  is  het  debiet,  dat  tussen  de  twee  platen  door  stroomt,  te  bereke- 
nen: 


AR 


t  .  2Lf 

h -  2  J  Lvrdr  =  -~ 

Jn  3  V 


(C.7) 


Deze  vergelijking  geeft  het  debiet  tussen  twee  platen.  Deze  relatie  kan  omgeschre- 
ven  worden  naar  de  stroming  tussen  kledingmateriaal  en  een  cilinder,  door  aan  te 
nemen,  dat  de  straal  van  de  cilinder,  R,  groot  is  ten  opzichte  van  de  luchtlaag,  2AR. 
In  dat  geval  geldt  namelijk: 

dx~Rd6  (C.8) 


Dit  geeft: 


.(«)=§£ 

3  TJ 


dPm 

(ARy 

de 

0; 

R 

(C.9) 


Het  debiet  verandert  als  functie  van  de  hoek,  omdat  lucht  het  kledingmateriaal  in  of 
uit  zal  stromen.  De  gemiddelde  snelheid  waarmee  de  lucht  tussen  het  kledingmate¬ 
riaal  en  de  cilinder  door  stroomt,  is: 


<  v(0)>= 


*,(?) 

2ARL 


1 

dPin 

\ 

3ij 

de 

0; 

R 

En  de  maximale  snelheid  (in  het  centrum  tussen  de  platen,  bij  r  =  0)  is: 


v(0)max 


1 

dP ln 

\ 

(ARf 

2jj 

de 

9) 

R 

(C.10) 


(C.11) 
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Bijlage  D  Analytische  oplossing  van  de  luchtsnelheid; 
laminaire  stroming 

Het  drukprofiel  onder  het  kledingmateriaal  volgt  uit  vergelijking  (D.l): 

2  L<  dX  V/M?)1  lOLflAP  10LR  (/,  p  ) 
dG2  gj  R  Rma,Pwg 

Indien  de  luchtstroom  om  de  kleding  laminair  is,  is  vergelijking  (D.l)  analytisch 
oplosbaar  (zie  ook  [5]).  In  bijlage  A  is  een  relatie  afgeleid  voor  de  buitendruk  als 
functie  van  de  aanstroomhoek  voor  lage  Reynolds-getallen  (Re  <  50).  In  dat  geval 
was  de  buitendruk  gelijk  aan: 

PA=P0  +^Pv'o(1-4 sin2  9)  (D-2 

Omdat  het  drukverschil  wordt  genomen  in  vergelijking  (D.l),  valt  Po  weg.  Als 
vervolgens  vergelijking  (D.2)  ingevuld  wordt  in  vergelijking  (D.l),  wordt  gevon- 
den: 

=  k{p uil  (0)(l  -  4 sin 2  d)-  Pin  )  (D.3 


30  77  ft2 

2  R^g  (ARj 


<D-5> 

Als  nu  aangenomen  wordt,  dat  de  binnendrak  een  functie  is  van  de  aanstroomhoek 
op  een  analoge  manier  als  de  buitendruk: 

Pin=A-B  sin20  (D-6) 

dan  kan  door  invulling  van  vergelijking  (D.6)  in  vergelijking  (D.3)  gevonden 
worden  dat: 


B=puiM- 
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D.2 


zodat: 
P 


(0)=  Pui,  (0)f  1 - - - sin2  #1 

v  ’  ""v  \  4  +  tf  4  +  A-  J 


(D.9) 


De  luchtsnelheid  (in  m/s)  door  het  kledingmateriaal  heen,  is  gerelateerd  aan  het 
drukverschil  (vergelijking  (D.10)): 


v  =  ■ 


10 


R matPwS 


(D.10) 


Beide  drukken  zijn  in  dit  laminaire  geval  bekend,  zodat  voor  de  luchtsnelheid  door 
het  kledingmateriaal  heen  gevonden  wordt: 


v  = 


80  PUM 
Rma,PwS  (4+^) 


(l-2sin20) 


(D.ll) 


Of  als  de  buitendruk  bij  hoek  0  ingevuld  wordt: 


v  =  Kh 


Waarbij  de  constante  Kb  gelijk  is  aan: 


Kb 


Pi  40 

Pwg( 4+^) 


(l-2sin2#) 


(D.12) 


(D.13) 


Bij  nadere  beschouwing  is  deze  constante  dus  niet  constant,  maar  een  fimctie  van 
de  aanstroomhoek.  Toch  is  vergelijking  (D.12)  op  deze  manier  weergegeven, 
omdat  hij  dan  goed  te  vergelijken  is  met  de  empirische  relatie  (D.14): 


v  =  K. 


(D.14) 


Relatie  (D.14)  geeft  de  gemiddelde  luchtsnelheid  door  het  kledingmateriaal  op  de 
plaatsen,  waar  de  luchtstroom  het  kledingmateriaal  in  gericht  is.  Om  relatie  (D.14) 
dus  theoretisch  te  verifieren,  moet  eerst  gemiddeld  worden.  Als  eerste  moet  geke- 
ken  worden,  tussen  welke  twee  hoeken  de  stroom  inwaartse  gericht  is.  De  empiri¬ 
sche  formule  ging  ervan  uit,  dat  de  stroom  inwaarts  gericht  is  tussen  0  =  ±45°.  Op 
de  grens  tussen  inwaartse  en  uitwaartse  luchtstroom  zal  gelden:  v=0.  Dus: 

l-2sin2  6  =  0  (D.15) 


Hieruit  volgt  inderdaad,  dat  0  =  ±45°. 

De  snelheid,  die  volgt  uit  vergelijking  (D.ll),  is  plaatsgebonden.  Om  een  verge- 
lijkbaar  getal  te  vinden  met  de  empirische  waarde  van  Ke,  zal  dus  de  snelheid 
gemiddeld  moeten  worden  tussen  -45°  en  45°.  In  formulevorm: 
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<  V> 


Pi  40  Vp 
8Pw  (4  +  ^)/?ma, 


4 


p,  40  Vp2 
£Pw  (4  +  K)  Rmat  [0f^ 

4 

(D.16) 


Dit  geeft  dus: 

Pi  40 


<  v>= 


gpw  (4  +  K)Rmal  n 


(D.17) 


Zodat  voor  de  constante,  Kb,  gevonden  wordt: 

K  80  Pi 
b  n{4+K)gpw 


(D.18) 


De  empirische  formule  (D.14)  werd  gevonden  door  gebruik  te  maken  van  de  vol- 
gende  experimentele  omstandigheden  [8]: 
vo  =9,1  m/s; 

2AR  -  2  *  10-3  m; 

R  =  25  *  lO-3  m; 

Rmat  =  2,66  mmH20/(cm7s). 

Dit  levert  op: 

Re  =  3  *  104; 

K  =  6  *  lO-3  en  Kb  =  7,66  *  10'4  (mmH20/(cm/s))/(m/s). 

Experimenteel  werd  gevonden:  Ke  =  5,7  *  10"4  (mmH20/(cm/s))/(m/s). 

De  laminaire  schatting  is  dus  circa  30%  te  hoog  ten  opzichte  van  de  experimentele 
waarde.  De  stroming  is  verre  van  laminair  in  dit  geval  (zie  Reynolds).  Een  verschil 
was  dus  te  verwachten. 


[5]  Fedele,  P.D., 

‘Model  of  aerosol  protection  offered  by  permeable  protective  garments’, 
American  Society  for  Testing  and  Materials,  Philadelphia,  (1992). 

[8]  Wal,  J.F.v.d., 

‘Permeabele  beschermende  kleding,  n  De  penetratie  van  aerosol  door  kle- 
ding’, 

TNO-rapport  1967-22,  (1967),  Rijswijk. 


TNO-rapport 


PML  1998-A107 
Bijlage  E 


E.1 


Bijlage  £  Analytische  oplossing  van  de  luchtsnelheid; 
turbulente  stroming 


Indien  de  luchtstroom  turbulent  is,  zal  de  drukverdeling  om  de  cilinder  anders  zijn 
dan  in  het  laminaire  geval.  Dit  heeft  dus  gevolgen  voor  de  snelheid  door  het  kle- 
dingmateriaal,  aangezien  aangenomen  wordt  dat  het  ontwikkelde  model  blijft 
gelden.  Om  plaatselijke  luchtsnelheden  te  berekenen  moet  uitgegaan  worden  van 
gemeten  drukprofielen.  In  het  turbulente  geval  is  dus  geen  analytische  oplossing  te 
vinden  voor  het  model.  Ge'integreerd  over  het  totale  oppervlak  van  de  cilinder  is 
wel  wat  bekend  over  de  buitendruk: 


Ku  —Cw  2  P'vo 


(E.1) 


Waarbij  Cw  een  wrijvingsterm  (‘drag’)  is.  Deze  term  is  een  functie  van  Reynolds  in 
het  laminaire  gebied,  maar  in  het  turbulente  gebied  wordt  hij  ongeveer  constant. 
Voor  cilinders  geldt:  Cw==  1.  Als  vervolgens  de  binnendruk  verwaarloosd  wordt  ten 
opzichte  van  de  buitendruk,  dan  kan  gevonden  worden: 


10  1  -  ,.a  _  5Pi  vo 

'  fJjVQ  -  —————— 

RnuuSPw  2  8PW  Rmn, 


(E.2) 


De  constante,  Kb,  is  dan  gelijk  aan: 

Kb  (E.3) 

8Pw 

Dit  geeft  Kb  =  6,0  *  10'4  (mmH20/(cm/s))/(m/s).  Dit  ligt  dus  erg  dicht  bij  de 
experimentele  waarde  van  Kg:  5,7  *  10'4  (mmH20/(cm/s))/(m/s). 

Een  nadeel  van  deze  rekenwijze,  is  dat  niets  gezegd  kan  worden  over  plaatselijke 
drukken  als  functie  van  de  aanstroomrichting.  Verder  kan  de  invloed  van  waarden, 
zoals  de  afstand  tussen  het  kledingmateriaal  en  de  cilinder,  op  de  luchtsnelheid  niet 
bepaald  worden. 
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Bijlage  F  Parameteranalyse 


F.l  Laminaire  stroming 


In  de  parameterstudie  is  steeds  per  parameter  gekeken  naar  het  effect  op  Kb-  Het 
zou  handig  zijn,  om  de  effecten  van  de  verschillende  parameters  samen  te  vegen  in 
een  nieuwe  parameter,  zodat  het  effect  van  alle  verschillende  parameters  vergele- 
ken  kan  worden  in  6en  plaatje.  In  het  laminaire  geval  is  wel  bekend  hoe  de  ver¬ 
schillende  parameters  invloed  uitoefenen  op  Kb-  Deze  relatie  wordt  afgeleid  in 
Bijlage  D  en  ziet  er  als  volgt  uit: 


7t\ 


4  + 


80 _ 

30  7]  ft2  " 

7  RmatPwS  (AR 


Pi 

8Pw 


Dit  is  ook  te  schrijven  als: 


±.=jlm+3„jl — *! — 

Kb  20  Pi  Pi  (2 PRfRm, 


(F.l) 


(F.2) 


Als  hierbij  de  bekende  waarden  ingevuld  worden,  dan  kan  het  volgende  rechtlijnige 
verband  gevonden  worden: 


—  =  1303/58+7216*  10~5 
Kb 


/g2 

(?pr 


(F.3) 


Dus  een  plot  van  de  reciproke  waarde  van  Kb  tegen  de  verhouding 
fg2 

(2 

geeft  in  het  laminaire  geval  een  rechte.  In  het  laminaire  geval  kan  de  in¬ 
vloed  van  de  verschillende  parameters  dus  samengevoegd  worden  tot  de¬ 
ze  verhouding. 


F.2  Turbulente  stroming 

Het  ligt  voor  de  hand  om  te  bekijken,  hoe  de  uit  het  drukprofiel  berekende  waarden 
zich  gedragen  als  functie  van  de  bovengenoemde  verhouding.  In  figuur  F.l  is  dit 
verband  te  zien. 
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RlA(2AR)*R«<)  ((anfcWranHjO  m)) 


Figuur  FJ:  De  reciproke  van  Kb  alsfunctie  van  een  gecombineerde  parameter. 

Uit  figuur  F.l  zijn  verschillende  resultaten  te  halen.  AIs  eerste  is  duidelijk  te  zien, 
dat  de  empirische  formule  (Van  der  Wal  in  figuur  F.l)  geen  rekening  houdt  met  de 
invloed  van  de  verschillende  parameters.  Er  is  een  constante  waarde  aangenomen. 
Ter  vergelijking  is  het  verband  weergegeven  voor  het  laminaire  geval  (vergelijking 
(F-3))- 

Verder  is  te  zien,  dat  deze  verhouding  een  goede  manier  is  om  een  totaalbeeld  van 
de  invloed  van  de  verschillende  parameters  te  krijgen.  Er  wordt  een  goed  verband 
gevonden  tussen  deze  verhouding  en  de  reciproke  van  Kb-  Hoewel  dit  verband  niet 
helemaal  lineair  is,  geeft  de  volgende  empirische  relatie  toch  een  correlatie- 
coefficient  van  r2  =  0,991: 


—  =  1722  +  10,85*  10~5 

Kb 


j?2 

{2ARyRmat 


(F.4) 


De  overeenkomst  tussen  vergelijking  (F.3)  en  (F.4)  is  opvallend.  Overigens  is  het 
sprekend  om  vergelijking  (F.4)  als  volgt  te  schrijven: 


Kh 


5,8*  10 ' 


1+63*10 


-8 


R A 


{2MyRMI 


(F.5) 


Als  de  invloed  van  de  verschillende  parameters  verwaarloosd  wordt  (het  tweede 
deel  van  de  noemer  van  vergelijking  (F.5)),  wordt  op  1 ,9%  nauwkeurig  de  empiri¬ 
sche  waarde  van  Ke  gevonden.  Dit  geeft  wel  erg  veel  ondersteuning  aan  zowel  de 
metingen  [8]  als  de  hier  uitgevoerde  theoretische  berekeningen.  Het  heeft  wel  als 
resultaat,  dat  de  empirische  formule  in  sommige  gevallen  iets  aangepast  moet 
worden.  Deze  aanpassing  kan  voor  lage  luchtweerstanden  (‘Out  Of  Area’- 


<■ 
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F.3 


materialen)  belangrijk  zijn.  In  dat  geval  zal  de  luchtsnelheid  door  kledingmateriaal 
berekend  moeten  worden  met: 


vmat 


5,8*10' 


vo 


1+63*10 


-8 


R A 


R 


mat 


(P.6) 


Zoals  al  eerder  opgemerkt  is,  geeft  deze  formule  alleen  een  waarde  voor  de  gemid- 
delde  snelheid  door  het  kledingmateriaal  op  de  plaatsen  waar  de  luchtstroom  naar 
binnen  gericht  is. 


[8]  Wal,  J.F.  van  der, 

‘Permeabele  beschermende  kleding,  II  De  penetratie  van  agrosol  door  kle- 
ding’, 

TNO-rapport  1967-22,  (1967),  Rijswijk. 
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